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做 流 变 胶 济 回 曲线 三 参数 模型 建 模 及 参数 识别 


陈 冬 
(河北 化 工 医药 职业 技术 学 院 ,050026 石家庄 ) 


摘 要 :利用 安 东 帕 流 变 仪 (MRC302 ) 测 试 了 兰 基 铁 粉 质量 分 数 为 40% 的 磁 流 变 胶 应 力 - 应 变 ,应 
力 -应 变速 率 两 个 灌 回 曲线 (频率 0.1 Hz, 应 变 幅 值 22% ,温度 25 ,电流 分 别 为 0 A、0.5 A、l A、 
2 A.3 A、4A,5 A)。 利 用 遗传 寻找 最 优 算法 ,建立 三 参数 液体 模型 ,对 磁 流 变 胶 的 应 力 -应 变 ,应 力 - 
应 变速 度 进行 参数 识别 。 Pa 与 实验 结果 的 均 方 差 确定 参数 模型 在 各 种 电流 下 对 磁 流 变 
胶 动 态 特 ， 性 建 模 的 满 满意 
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Numerical simulation and analysis of transmission characteristics 
of safety coupling based on MR gel under Halbach 
permanent magnet arrangement 


CHEN Dong 


( Hebei Chemical & Pharmaceutical Vocational Technology College ,050026 Shijiazhuang, China) 


Abstract :In this paper, hysteresis curves (frequency ,0. 1 Hz;strain amplitude ,22% ;temperature ,25 C; 
current,0 A,0.5A,1 A,2A,3A,4A,5 A,respectively) of stress-displacement and stress-velocity of mag- 
netorheological (MR) gel were measured using MRC302. A three-parameter liquid model was established 
to identify the parameters of stress displacement and stress velocity of MRF by using genetic algorithm. Ac- 
cording to the mean square deviation of the fitting results and the experimental results ,the satisfaction of the 
parameter model to the dynamic characteristic modeling of MR gel under various currents was determined. 
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磁 流 变 胶 是 一 种 新 型 磁 流 变 智能 材料 ,是 由 软 ” 胶体 状 的 流动 状态 ,一 旦 施加 磁场 ,会 在 毫秒 级 时 间 


磁 微米 级 (3 ~5_ hm) 颗粒 (如 铁 粉 、 普 基 铁 粉 等 ) ， 
均匀 悬浮 在 非 导 磁 基 液 “ (如 聚氨酯 、 硅 树脂 
等 ) 中 形成 的 一 种 智能 复合 材料 。 磁 流 变 胶 的 流 变 


特性 与 磁 流 变 液 相似 , 即 在 无 磁场 环境 下 ,其 表现 出 
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内 了 迅速 变 成 类 固态 ,并 且 具 有 一 定 的 剪 切 届 服 应 力 
(和 剪 切 屈服 应 力 随 着 磁场 强度 增 大 而 先 增 大 后 趋 于 
一 个 饱和 值 ) , 当 撤 下 磁场 后 又 迅速 恢复 原来 的 胶体 
状态 。 整 个 过 程 具有 快速 .可 逆 \ 可 控 的 特点 ”。 磁 


力学 学 报 ,2022 ,39(5) :932-937. 
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流 变 胶 由 于 其 基 液 是 黏度 比 硅油 更 高 的 聚氨酯 或 者 
硅 树 脂 等 高 分 子 材料 ,因此 其 基体 黏度 比 磁 流 变 液 
更 大 ,能 够 防止 颗粒 由 于 久 置 而 发 生 沉降 现象 "8 。 
因此 , 磁 流 变 胶 在 某 些 工程 应 用 中 能 够 代替 磁 流 变 
液 发 挥 更 大 的 优势 ,如 磁 流 变 阻 尼 器 "1 、 磁 流 变 离合 
器 59 \ 磁 流 变 吸 振 器 和 磁 流 变 制动器 等 。 

磁 流 变 胶 的 动态 力学 性 能 预计 介 于 磁 流 变 液 和 
磁 流 变 弹 性 体 之 间 ,在 施加 激励 时 表现 出 黏 性 和 /或 
弹性 行为 。 正 如 预期 的 那样 ,对 于 磁 流 变 材料 动态 
特性 ,建立 计算 效率 高 的 模型 ( 即 参数 较 少 的 模型 ) 
准确 地 描述 动态 力学 响应 是 其 工程 应 用 的 基础 。 文 
献 [13] 提 出 了 秋 - 弹 -塑性 模型 来 预测 MR 阻尼 器 的 
动态 消 回 特性 。Bouc-Wen 注 回 算 子 用 于 描述 材料 
或 器 件 中 力 -变形 关系 的 固有 灌 回 特性 。 为 了 更 
准确 地 描述 磁 流 变 阻 尼 器 在 大 范围 输入 (振幅 .频率 
E ee 


uc-Wen 模型 。 考 虑 到 磁 流 变 材料 在 纯 剪 切 模 
式 环 的 剪 切 确 化 和 惯性 效应 ,文献 [16] 提出 了 基于 
Boie-Wen 清 回 模型 的 无 弹 秘 单元 动力 学 模型 。 

一 文献 [17] 提 出 准 静 态 一 维 模型 ,该 模型 假设 磁 
框 硬 作用 仅仅 存在 于 链 中 2 个 相连 颗粒 间 ; 并 且 应 
i 
设 5 将 模型 运用 到 磁 流 变 弹 性 体 中 得 到 较 好 的 拟 合 
结 洪 。 文 献 [18] 建 立 了 考虑 链 中 所 有 颗粒 影响 的 

模型 ,该 模型 根据 磁 弹 性 力学 的 基本 最 小 能 量 
原理 ,计算 磁化 强度 分 布 并 对 磁 流 变 弹性 体 的 力学 
性 老 进 行 拟 合 。 文 献 [ 19] 针对 磁 流 变 液 的 力学 特 
性 提出 了 双 笑 度 模型 ,该 模型 根据 外 部 剪 切 应 力 的 
大 名 将 力学 行为 分 为 届 服 前 阶段 和 届 服 后 阶段 ,并 
将 该 模型 运用 到 磁 流 变 液 中 ,结构 显示 该 模型 更 加 
适合 磁 流 变 液 工作 在 挤 压 模式 下 的 流动 情况 。 文 献 
[20] 在 Binham 模型 的 基础 上 提出 Herschel-Bulkley 
模型 ,该 模型 能 够 将 届 服 后 的 磁 流 变 液 的 剪 切 稀 化 
表现 出 来 ,能 够 比 Bingham 模型 更 加 准确 地 描述 磁 
流 变 液 的 力学 模型 。 

以 上 模型 都 是 针对 磁 流 变 胶 以 外 的 材料 ( 磁 流 
变 液 和 磁 流 变 弹 性 体 ) 进行 建 模 和 参数 识别 ,而 磁 流 
变 胶 作为 新 型 磁 流 变 智能 材料 ,在 工程 应 用 中 具有 
潜在 应 用 价值 。 因 此 ,针对 磁 流 变 胶 力学 特性 的 建 
模 研究 具有 一 定 的 价值 和 意义 。 基 于 磁 流 变 胶 离合 
器 的 输出 扭矩 可 控 , 其 优点 在 于 能 够 实现 半 主 动 控 
制 。 而 实现 控制 的 前 提 就 是 需要 掌握 磁 流 变 材料 的 
动态 流 变 特 性 ,并 对 其 在 各 个 激励 条 件 下 的 流 变 响 
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应 展开 参数 化 建 模 , 建立 一 般 化 的 能 够 准确 有 效 预 
测 磁 流 变 胶 动 态 流 变 行为 的 参数 模型 ,为 实现 磁 流 
变 絮 件 半 主 动 控制 打下 理论 基础 。 

本 研究 主要 针对 磁 流 变 胶 的 应 力 - 剪 切 应 变 、 应 
力 - 剪 切 速率 的 清 回 曲线 , 基于 遗传 优化 算法 ,建立 
三 参数 液体 模型 ,并 对 其 展开 力学 特性 的 参数 建 模 
和 参数 识别 。 最 后 根据 拟 合 结果 与 实验 结果 的 均 方 
差 , 对 参数 模型 在 不 同 电流 下 对 磁 流 变 胶 动态 特性 
建 模 的 满意 度 。 


1 磁 流 变 胶 滞 回 曲线 测试 


测试 装置 如 图 1 所 示 。 研 制 狭 基 铁 粉 质量 分 数 
为 40% 的 聚氨酯 基 磁 流 变 胶体 如 图 1(a) 所 示 。 在 
室温 条 件 (25 C ) 下 对 研制 的 磁 流 变 胶 应 力 -应 变 、 
应 力 -应 变速 率 曲 线 进行 测试 。 图 1(a) 所 示 测 试 工 
具 为 安 东 帕 (型 号 :MRC302 ) 流 变 仪 。 流 变 仪 两 测 
试 平板 直径 均 为 20 mm, 设置 流 变 仪 两 平板 之 间 的 
距离 为 1 mm。 图 1(b) 和 图 1(c) 分 别 为 磁 流 变 胶 在 
无 外 加 磁场 下 (表现 出 类 胶体 ) 和 外 加 磁场 下 (表现 
出 类 固体 特性 ) 的 宏观 现象 。 测 试 条 件 为 : 频率 
0.1 Hz ,应 变 幅 值 22% ,施加 电流 0 A、0.5 A、l A、 
2A3A4A.5A, 


IN 


(c) 施加 磁场 


聚氨酯 磁 流 变 胶 
(a) 


妈 1 测试 装置 
Fig.1 Test device 
不 同 电流 下 应 力 -应 变 、 应 力 -应 变速 率 如 图 2 所 
示 。 从 图 2(a) 中 可 以 看 出 , 当 施 加 电流 为 0 A 时 ， 
磁 流 变 胶 应 力 -应 变 表现 出 黏 性 特性 , 当 施加 电流 
后 , 磁 流 变 胶 应 力 - 应 变 曲线 表现 出 具有 黏 弹性 特 
性 ,并 且 黏度 ( 滞 回 曲线 包 络 的 面积 ) 和 弹性 ( 滞 回 
曲线 最 大 应 力 和 最 小 应 力 之 间 的 连接 线 ) 随 着 外 部 
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电流 强度 增 大 而 增 大 。 这 是 由 于 随 着 电流 的 增 大 ， 
磁 流 变 胶 所 在 区 域 磁场 强度 越 大 ,成 链 越 粗 导致 的 。 
从 图 2(b) 中 可 以 看 出 ,施加 电流 后 磁 流 变 胺 的 沾 回 
曲线 具有 明显 的 非 线性 特性 。 
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> 图 2 磁 流 变 胶 应 力 -应 变 应 力 -应 变速 率 曲线 


NE 2 Stress-strain and stress-strain rate curves of MR gel 


®& 
2 三 三 参数 模型 
三 
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由 于 磁 流 变 胶 在 无 磁场 下 为 一 类 黏 性 流体 ， 
此 其 动力 学 模型 需要 一 个 初始 的 黏 性 单元 ci ,在 施 
加 磁场 后 ,表现 为 类 固体 状态 ,具有 一 定 的 磁 控 黏 弹 
性 特性 , 即 当 简 谐 应 变 作用 在 其 应 变 -应 力 沾 回 曲线 
时 ,同时 具有 黏 性 特性 和 弹性 特性 。 因 此 ,在 磁场 作 
用 下 ,其 动力 学 模型 需要 用 可 控 的 黏 弹性 单元 表示 ， 
这 部 分 采用 Kelvin 模型 单元 中 由 黏 性 单元 c, 和 弹性 
单元 表示。 男 外 ,同时 需 能 够 准确 表示 其 在 有 无 
人 磁场 下 的 忒 弹性 行为 ,这 一 部 分 用 c, 来 表示 。 因 此 ， 
其 总 体 参数 模型 如 图 3 所 示 ( 称 为 三 参数 模型 ) 。 

图 3 中 三 参数 液体 模型 中 有 Kelvin 模型 单元 
(图 中 虚线 部 分 ) 。Kelvin 模型 单元 中 有 黏 性 单元 c 
和 弹性 单元 h。 从 图 3 中 可 以 看 出 三 参数 液体 模型 
是 在 Kelvin 模型 单元 串联 一 个 黏 性 单元 c 。 因 此 ， 
三 参数 液体 模型 可 以 通过 式 (1) 描 述 , 即 
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下 二 万 =cixi 
| (1) 


F+f=hk(x—-x) +c (x—Xi) 
式 中 :有 为 模型 应 力 ;及 为 初始 应 力 。 


Kelvin 模型 
图 3 三 参数 液体 模型 
Fig.3 Three-parameter fluid model 


模型 计算 所 得 应 力 与 实测 应 力 之 差 的 均 方 根 值 
与 对 应 电流 下 实测 的 最 大 应 力 比 值 作 为 模型 满意 
度 ,如 式 (2) 所 示 。 


1 = 7 x 100% (2) 


式 中 :作为 模型 满意 度 ; V 为 数据 个 数 ;下 ,为 第 n 个 
所 的 模型 预测 值 ;7, 为 第 n 个 点 的 实测 值 ;WW 为 各 个 
电流 下 实测 的 最 大 应 力 值 。 


3 结果 与 讨论 


在 不 同 电流 激励 下 ,利用 遗传 算法 对 该 模型 参 
数 进 行 最 优 识别 ,具体 优化 算法 流程 如 图 4 所 示 。 


设置 模型 参数 和 优化 
边界 条 件 

初始 化 种 群 

个 体 适应 度 评价 


痊 出 最 优 模型 参数 


色 体 的 群 


图 4 GA 最 优 求解 模型 参数 流程 示意 图 


Fig.4 Flow chart of GA optimal solution model parameters 


图 5 和 图 6 分 别 表示 不 同 电流 下 磁 流 变 胺 应 
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第 5 期 
力 -应 变 ,应 力 -应 变速 率 三 参数 液体 模型 拟 合 和 实 
测 值 对 比 曲线 。 从 图 5 和 图 6 可 以 看 出 ,三 参数 液 
体 模型 能 够 很 好 地 拟 合 磁 流 变 胶 应 力 -应 变 、 应 力 - 
应 变速 率 曲线 。 不 同 电流 下 三 参数 液体 模型 参数 识 
别 结果 如 表 1 所 示 。 
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到 5 不 同 电 流下 人 磁 流 变 胶 应 力 -应 变 

三 参数 液体 模型 拟 合 和 实测 值 对 比 曲 线 


Fig.3 Stress-strain comparison curves of three-parameter 


liquid model fitting and measured value of 


MR gel under different currents 


从 表 1 可 以 看 出 ,c, 和 c 的 值 随 着 电流 强度 增 
大 而 增 大 ,说 明 材 料 的 黏 性 随 着 电流 增 大 而 增 大 ,该 
结论 与 图 2(a) 中 的 清 回 曲线 面积 随 着 电流 增 大 而 
增 大 相符 。k 的 绝对 值 随 着 电流 强度 增 大 而 增 大 ， 
该 结论 也 与 图 2(a) 中 的 灌 回 曲线 平均 斜率 ( 灌 回 曲 
线 最 大 应 力 和 最 小 应 力 之 间 的 连接 线 ) 随 着 外 部 电 
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多 6 不 同 电流 下 磁 流 变 胶 应 力 -应 变速 率 
三 参数 液体 模型 拟 合 和 实测 值 对 比 曲 线 


Fig.6 Stress-rate comparison curves of three-parameter 


liquid model fitting and measured value 


of MR gel under different currents 


流 强度 增 大 而 增 大 相符 。 在 不 同 电 流下 , 三 参数 液 
体 模型 对 实测 值 的 拟 合 满意 度 均 小 于 10% ,在 接受 
范围 之 内 。 

三 参数 模型 中 的 3 个 参数 cc, 和 都 随 着 外 
加 电流 表现 出 单调 变化 ;具体 来 说 ,c, 和 c, 随 着 电流 
增强 而 单调 增 大 , k, 随 着 电流 增强 而 单调 减 小 。 


和 
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这 种 模型 参数 值 随 电流 单调 变化 有 助 于 后 期 控制 系 

统 的 设计 ,说 明 用 该 模型 描述 磁 流 变 胶 动 态 流 变 特 性 

具有 较 好 的 可 控 特 性 。 由 于 磁 流 变 材料 及 其 功能 


件 ( 如 阻尼 器 ) 具 有 内 在 非 线性 滞 回 特性 ,如 何 建 立 适 
ee 目前 面临 的 问题 


0 E 耗 型 屁 件 ,需要 对 材料 的 滞 回 曲线 


,应 中) 的 化 情况 ,建立 与 其 


研究 分 析 , 把 握 其 茸 弹 性 特性 随 外 部 激励 (如 电 


tt 浪 则 轩 信 eh 


。 因 此 ,三 参数 液体 模型 可 作为 磁 流 变 胶 特 性 的 ”应 的 参数 模型 进行 描述 ,具有 一 定 的 理论 指 
表 1 模型 参数 识别 结果 与 对 应 的 满意 度 
Tab.1 Model parameter identification results and corresponding satisfaction 
模型 电流 /A a c 局 hn 1H/ % 
三 参数 液体 模型 0 4.46 3.79 -6.33 13.46 7.33 
0.5 51.75 41.79 —9.84 282.60 8.11 
1 137. 15 111.09 —23.82 414.61 6.67 
2 280. 04 215.9%2 一 69.39 342. 64 6.93 
3 378.67 320. 49 -173.11 361.32 6.64 
4 470. 95 385.51 —305. 82 240.30 6.77 
5 567.27 429.98 —320.93 18.21 7.49 
cillation shear tests[ J ]. Royal society of chemistry,2019 ,18(9 ) : 
10124-10134. 
[3] MEHARTHAJ H,SIVAKUMAR S M,AROCKIARAJAN A. Signifi- 
和 cance of particle size on the improved performance of magnetorheo- 
分 x 
本 研究 对 自 和 行 人 饶 制 的 狼 基 铁 粉 质 量 分 数 ， logical gels [J ]. Journal of magnetism and magnetic materials, 
40% 的 聚氨酯 基 磁 流 变 胶 的 滞 回 特性 进行 建 模 。 2019 ,490(15) : 165483. 
验 测 导 磁 流 变 胶 应 力 -应 变 ,应 力 - 应 变速 率 2 个 灌 [4] DE SOUSA SR G,DOS SANTOS M P,BOMBARD A J F. Magne- 
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下 渍 | 下 下 下 3 D 
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明 ,不同 电流 下 利用 三 参数 液体 模型 拟 合 后 的 满意 nad 和 et oscillatory shear[ J]. Journal of 
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私 均 小 于 10% 三 参数 型 能 够 4 
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